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Using a combination of local spin density and Hubbard 1 approximations we study the mechansim 
of exchange interacion in EuX (X=0, S, Se and Te). We reproduce known experimental results about 
bulk modulus, critical pressure for structural phase transition, magnetic ordering temperature, spin- 
wave dispersions as well as momentum- and tempearuture-dependent band shift. Our numerical 
results show pressure induced competition between the hybirization enhanced exchange interaction 
and Kondo-like coupling in EuO. Possible ways to enhance T c are discussed. 
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Europium monochalcogenides EuX (X=0, S, Se and 
Te) have been studied extensively since 1961 [l| . As the 
only known examples of Heisenberg ferromagnetism in 
nature, EuO and EuS have their Curie temperatures (T c ) 
of 69.15 K and 16.57 K, respectively. On the other hand, 
EuSe has a complex magnetic structure at low temper- 
atures, and EuTe is antiferromagneticjlj. It was found 
that doping of EuO by electrons results in 100% spin po- 
larization of the conduction electrons Q, and the material 
has a colossal magnetoresistance effect stronger than fa- 
mous manganites. Moreover, very recently EuO has been 
integrated with Si and GaN, making it very attractive for 
spintronic applications H, and the interest to these sys- 
tems has been renewed [5Tjl2|. 

As finding ways to raise the T c up to room temper- 
atures in EuX is of both fundamental and technologi- 
cal importance, many past studies of their magnetic ex- 
change mechanism appeared in the literature. Based on a 
model calculation, Kasuya [l3| proposed that the nearest 
neighbor exchange coupling Ji is induced by the indirect 
exchange between 4/ and 5d electrons of Eu while the 
superexchange between the 4/ states of Eu and p elec- 
trons of anion can be ignored. Liu and Lee [13] claimed 
based on their band structure calculation that anion va- 
lence band has an important contribution to both J\ 
and Ji. On the other hand, based on density func- 
tional calculation and Wannier function analysis, Kunes 
et.al 15[ emphasized the importance of hybridization be- 
tween the 4/ of Eu and 2p of O and the associated su- 
perexchange interaction. There have also been other the- 
oretical models [l6| to describe magnetic exchange mech- 
anism in EuX including the s-f model[l7|. Experimen- 
tally, optical spectroscopy finds a considerable 4/-5<i mix- 
ing and suggests the importance of f-d exchange [18j; the 
Mossbauer experiment emphasizes the effect of 6s band 
of Eu and supports the s-f model [lj|; the neutron diffrac- 
tion stresses the contribution from the anion p shells [20j; 
the angle resolved photoemission (ARPES) experiment 
observes a momentum-dependent temperature-induced 
band shift [8[ and contributions from Eu 4/— 5d exchange 
to Ji and from Eu 4/-0 2p exchange to J2 while the x- 
ray absorption spectroscopy Q claims that the exchange 



interaction is due to f-d mixing without involvement of 
the anion p states. It was also found that pressur e [jj , 
[25j], epitaxial strain (2|| [27| and carrier doping [il [28| can 
vary the T c of EuX significantly although still did not 
reach the values comparable with the room temperature. 

In this work we address the controversial issue of un- 
derstanding magnetic exchange mechanism in europium 
monochalcogenites using a recently developed linear re- 
sponse approach [29J which is based on a combination of 
density functional and dynamical mean field theories [3(| ■ 
We reproduce major experimental results regarding their 
spin wave dis per sions [31j , pressure dependent transition 
temperatures |20l |24| , and temperature dependent band 
structures [8[ . We provide conclusive theoretical insights 
to various contributions to magnetic exchange interac- 
tions. 

Our electronic structure calculations with the full po- 
tential linearized-muffin-tin-orbital (LMTO) method |32| 
are done using the local spin density approximation 
(LSDA) for the conduction electrons and atomic Hub- 
bard 1 (Hubl) self-energy to approximate localized na- 
ture of the Eu f-electronsjlO] with the on-site Coulomb 
interaction parameters U =7 cV and J =1.2 eV[6, 7]. 
We also check that our results are robust within the 
reasonable range of U's from 6 to 9 eV. With the elec- 
tronic structure information, one can evaluate the mag- 
netic interaction J in the Heisenberg Hamiltonian H = 
JijSi ■ Sj, based on a magnetic force theorem [33| 
that evaluates linear response due to rotations of mag- 
netic moments [29}. This technique has been used success- 
fully for evaluating magnetic interactions in a series of 
Mott insulating oxides p9|. cuprate[34| and pnictide[35[ 
superconductors . 

The ground state properties predicted by our 
LSDA+Hubl calculation including magnetic moments 
and energy gaps are found to be in agreement with ex- 
periment. Moreover, the obtained exchange splitting of 
conduction band is about 0.65 eV, which is close to the 
experimental value 0.60 eV Q. Since at ambient pressure 
EuX compounds crystallize in rock-salt structure, but 
change to CsCl-type structure at high pressures [22| . [23j . 
we perform our calculations for both NaCl and CsCl type 



2 



TABLE I: Comparison between calculated using 
LSDA+Hubl method and experimental values for the 
bulk modulus Bo, its first derivative B'q and phase tran- 
sition pressure P c in Europium monochalcogenides. The 
experimental values are given in parentheses. 



EuO 


EuS 


EuSc 


EuTe 


Bo(GPa) 105 (118 a ) 


61 (61 6 ) 


53 (52 b ) 


43 (40 b ) 


B' 3.2 (2.2 a ) 


2.8 


2.8 


2.8 


P c (GPa) 48 (47 a ) 


26 (22 6 ) 


17 (15 b ) 


14 (ll 6 ) 



structures for a number of different volumes. The to- 
tal energy vs volume, E(V), curves were fitted by the 
Murnaghan equation of state (EOS), and the obtained 
bulk modulus (Bo) together with its first derivative (B' Q ) 
are listed in Table I. The crystal phase stability is an- 
alyzed by evaluating the enthalpy (H = E + PV), and 
the phase-transition pressure P c is evaluated from cross- 
ing the H(P) curves. It is found that LSDA alone sig- 
nificantly overestimates the values of Bq and P c . After 
inclusion of the correlation effects, our LSDA+Hubl cal- 
culation reproduces the experimental values successfully. 
This is shown in Table I. The parameters of EOS and P c 
are not sensitive to IPs in the range 6—9 eV. 

Using our LSDA+Hubl method and the magnetic 
force theorem, we subsequently evaluate the exchange 
constants as the integral over the q space using (8,8,8) 
grid. The magnetic interactions in these materials are 
usually characterized by two exchange constants, J\ and 
J 2 , although there is an argument about longer range 
interactions 36] . Our linear response approach allows us 
to evaluate J's in real space, and the numerical results 
confirm that they are short range with the magnetic cou- 
pling further than the second nearest neighbor to be al- 
most equal to zero. As shown in Table II, our J\ and 
J2 are smaller than the values extracted from neutron 

while quite close to the ther- 
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scattering experiment 
modynamic data37|. Consistent with the experiment, 
we obtain that moving from EuO to EuTe, the strength 
of J 1 decreases, meanwhile J 2 changes from FM to AFM- 
like. For EuSe, the magnitude of J\ is smaller than J2, 
but notice that for the rock-salt structure the number of 
nearest neighbors and of second nearest neighbors is 12 
and 6, respectively, so there is a competition resulting in a 
complex magnetic structure at low temperatures 0- For 
EuTe, both J\ and J2 become AFM-like and the ground 
state chang es to AFM0. Using the mean-field approx- 
imation (MFA), we also estimate the magnetic ordering 
temperature. The agreement between the experimental 
and calculated T c is good, with the expected overestimate 
of the mean-field value. 

Based on the obtained exchange interactions we eval- 
uate the spin-wave dispersions of EuO along major high 
symmetry directions. This is shown in Fig.l where for 
comparison we also plot by symbols the results of neu- 
tron scattering measurements [3l| . Our numerical data 



TABLE II: Calculated and experimental nearest neighbor, 
Ji, and second nearest neighbor, J 2 , exchange couplings as 
well as magnetic transition temperature for EuX (X=0, S, 
Se and Te) in units of K. The positive/negative signs denote 
ferro/antiferro magnetic coupling and Curie/Neel tempera- 
ture. 

EuO EuS EuSe EuTe 



Our results J L 0.60 0.12 0.10 

J 2 0.03 -0.10 -0.18 

Thermodynamic [37] Ji 0.67 0.19 0.13 

J 2 -0.06 -0.08 -0.12 

Neutron Scattering[31] Ji 0.61 0.24 

J 2 0.12 -0.12 

Our results T c 81.1 19.6 -5.9 

Experimental data [2] T c 69.3 16.6 -7.1 
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FIG. 1: Calculated spin-wave dispersions of EuO in compar- 
ison with the experiment (circles) [3lJ] . 



are found in good agreement with the experiment which 
is done on polycrystalline samples thus accessing direc- 
tion averaged spin- wave excitations. 

One of the controversy about the mechanism of ex- 
change interaction is the superexchange via the p or- 
bital of anion which is believed to be negligible due to 
its small hybridization with localized 4/ orbit al[i~3l [T3|. 
However this viewpoint has been argued due to recent 
theoretical works [3, SB] and a very recent observation of 
the momentum-dependent shift of 4/ states of Eu and 
2p states of 0@|. It has even been suggested that the 
f—p superexchange is the leading factor which induces 
change of sign in J 2 when moving from EuO to EuTe[26|. 
To clarify this effect, we performed the calculation with 
turned off spin polarization of the conduction band by 
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FIG. 2: Temperature dependent LSDA+Hubl band structure 
of EuO. (a) Temperature is 5 K. The red and black lines 
denote the majority and minority spin state respectively, (b) 
temperature is above T c . By arrows we show the experimental 
spin splitting in the 2p O band and the momentum dependent 
shift for the top of the 4/ Eu band. 



using LDA instead of LSDA approach. It is interesting 
that such LDA+Hubl calculation still gives a consider- 
able spin splitting for the p band of anion, and the mag- 
nitude of this splitting is almost the same as we obtain 
with the LSDA+Hubl method. Within the LDA+Hubl 
framework this spin splitting of the p band can only come 
from the hybridization with the 4/ band of Eu. Thus, we 
confirm that the f-p overlap is not negligible and the 
spin polarization in 4/ will result in shifting p band as 
suggested by the recent experimental work[8j. 

To further understand the effect of the f-p hybridiza- 
tion and in order to make comparisons with recent 
ARPES data Q regarding the temperature and moment- 
dependent Eu-4/ and 0-2p band shifts we perform the 
LSDA+Hubl calculation at temperatures below and 
above T c . Agree with the experiment 0, our calcula- 
tion shows that at low temperature (5 K), 0-2p state 
have considerable spin splitting as shown in Fig. 2 (a). 
With increasing temperature the spin splitting decreases 
and eventually becomes zero when temperature becomes 
above T c as shown in Fig2.(b). We also reproduce the 
temperature and moment-dependent Eu-4/ band shift. 
As shown in Fig. 2, our calculation predict that from low 
temperature to high temperature, the top of the Eu-4/ 
band at T point and X point have the shift of 0.35 and 
0.18 eV, respectively. These values are in good agreement 
with the experimental values of 0.32 and 0.07 eV[8|. 

Based on these arguments and noticing that EuX has 
the NaCl-type structure with anion located between two 
Eu atoms, one can then expect that the 4/-2p-4/ su- 
perexchange should play an important role in the ex- 
change mechanism. However, our values of J\ and Ji 
extracted from LDA+Hubl calculation are coming out 
to be very small (e.g., for EuO, both J\ and J2 < 0.001 



K) as compared to LSDA+Hubl result shown in Table 
II. This emphasizes that the major part of the exchange 
process is going thru virtual excitations to the conduc- 
tion band. If one does the same test for Mott insulating 
oxides[29| and cuprates[34j where the 3d-2p-3d superex- 
change dominates, one finds that the results of the cal- 
culations with LDA+Hubl and LSDA+Hubl are indeed 
similar to each other. This clearly indicates that the spin 
polarization in the conduction band is essential to explain 
the magnetic behavior of EuX. At the same time one can 
also conclude that there is considerable f-p hybridiza- 
tion in EuX system as suggested by recent experiment @ 
and theory work [lj| (26[. However the f-p superexchange 
being the second order effect can be ignored for those sys- 
tems. 

To see which exactly orbitals are contributing to the 
exchange process, we additionally made a calculation of 
J's with an artificial external potential applied to a par- 
ticular orbital. It turns out that a downshift of bp or- 
bital for Eu or of s orbitals for anion does not affect the 
exchange constants. On the other hand, they are ex- 
tremely sensitive to the position of 5 c? states of Eu: even 
a small up-shift here significantly decreases the strength 
of Ji and J2 due to the increase of band separation be- 
tween 4/ and 5d and following 5d-4f dehybridization. 
Shifting 6s level of Eu affects the J\ and J2 as well, al- 
though the effect is smaller than for 5d. Since anion 
atoms are between second nearest neighbor Eu ions, it is 
natural to expect that J2 is mediated by the 2p electrons 
of anion[§|, |20| through the d-p and s-p hybridization. 
As a result the value of J2 should be sensitive to the po- 
sition of the 2p band. Our numerical calculation, on the 
other hand, shows that both J\ and J2 are almost insen- 
sitive to the shift of the p band of anion. This result is 
consistent with our previous conclusion on the smallness 
of 4/-2p-4/ superexchange. For example, we upshift the 
p level of O in EuO by 1.0 eV but J\ and J 2 are almost 
unchanged. So, our calculation clearly rules out possible 
contributions from the d-p and s-p hybridization, and 
even the second nearest neighbor coupling J2 is mediated 
by the 6s and 5d electrons. 

We have studied the effect of pressure which does not 
only enhance the hybridization between the conduction 
band and the / states, but also increases the crystal-field 
splitting between 5d t2 g and e g states of Eu and reduces 
the energy gap between the t 2g and the 4/ states [20|. 
Consistent with the experiment, our theoretical results 
show that pressure enhances the exchange constants and 
T c of EuS, while it changes the ground state of EuSe and 
EuTe from AFM to FM. However, for EuO our numerical 
data agree with the experiment only at the low pressure 
region, while for pressures larger than 20 GPa, our calcu- 
lated Curie temperature still rises but the experimental 
one decreases [2J. This is illustrated in Fig. 3. Our calcu- 
lation shows that pressure closes the band gap and results 
in metallic behavior of EuO, but even for pressures as 
high as 40 GPa, the calculated exchange coupling is still 
short range. Based on our LSDA+Hubl calculation, we, 
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FIG. 3: Pressure dependence of magnetic transition temper- 
ature of EuO experimental (triangle) [3] , theoretical(circles) 
as well as Pressure dependence of Kondo coupling Jk- 



on the other hand, estimate the on-site Kondo cou pling 
strength Jk using a method described by us earlier |38j. 
As shown in Fig. 3, Jk increases rapidly above 20 GPa, 
which results in not only the Kondo-like screening but 
also in an AFM like intersite coupling in the second order 
perturbation theory with respect to Jk- Both effects will 
suppress the T c as is seen experimentally [24|. 

We have finally studied the effect of electron doping 
on T c for EuO using the virtual crystal approximation 



(VCA). Our calculation shows that there are two com- 
peting factors: since the bottom of the conduction band 
consists mainly of the majority spin, the doped electron 
will enter the spin-polarized manifold, and this results in 
the onset of moment in 5d band which will enhance T c . 
On the other hand, free carriers will induce the RKKY 
interaction which will supress T c . We indeed find that 
our theoretical T c first increases as the function of dop- 
ing , and then goes thr oug h a maximum in accord with 
the experimental trend[28(. So we conclude that a com- 
bination of both doping and pressure can be the efficient 
ways to reach higher T c in those systems. 

In summary, we have calculated the exchange con- 
stants of EuX and reveal contributions from various or- 
bitals relevant to the exchange mechanism. In particular, 
we showed that the 5d-Af and 6s-4/ indirect exchange 
is dominant while the p electrons of anion has no con- 
tribution to it. We also reproduced that bulk modulus, 
pressure-induced phase transition and magnetic ordering 
temperatures as well as the spin-wave dispersions agree 
well with the experiment. We also suggest that the ab- 
normal behavior of pressure dependent T c in EuO is due 
to increase in Kondo-like coupling in the metallic phase. 
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